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Abstract 

The  motivation  of  this  work  is  quantify  the  degra¬ 
dation  of  aging  electrical  cables.  The  dielectric  ma¬ 
terial  parameter  of  insulation  can  be  correlated  with 
degradation.  In  this  paper,  the  forward  problem  is 
posed  as  a  microwave  nondestructive  evaluation  (NDE) 
problem.  A  2D  finite  difference  time  domain  (FDTD) 
was  used  to  evaluate  the  forward  problem.  The  in¬ 
verse  problem  is  solved  via  a  Bayesian  approach.  The 
Bayesian  formulation  describes  the  solution  as  a  pos¬ 
terior  distribution  over  the  model  parameters  such  as 
complex  permittivity.  Since  there  is  no  analytical  so¬ 
lution  for  the  posterior  distribution,  the  Markov  Chain 
Monte  Carlo  (MCMC)  method  is  employed  to  numer¬ 
ically  solve  for  it.  The  Metropolis-Hasting  algorithm 
is  used  in  particular.  Results  for  computational  ex¬ 
periments  are  demonstrated  to  show  feasibility  of  this 
approach. 

1  Introduction 

The  degradation  of  insulation  in  electrical  cables  used 
with  instruments  and  control  facilities  is  of  serious  con¬ 
cern  for  management  of  aging  large  scale  systems  such 
as  airplanes,  power  plants,  etc  [1].  The  fact  that  di¬ 
electric  properties  of  materials  used  in  electrical  cables 
change  as  they  degrade  over  the  course  of  their  ser¬ 
vice  life  provides  a  means  of  quantifying  the  degrada¬ 
tion  if  the  dielectric  property  can  be  measured  accu¬ 
rately.  Microwave  NDE  can  make  these  measurements 
and  thus  quantitatively  assess  the  degradation  of  cable 
insulation.  The  idea  has  been  already  applied  to  the 
in-service  inspection  in  a  variety  of  industrial  applica¬ 
tions.  Time-Resolved  Microwave  Dielectric  Absorption 
(TRMDA)  enables  the  measurement  of  the  evanescent 
electrical  field  leaking  from  the  probe  into  the  cable 
insulation  by  using  cavity  resonator  [2].  This  electrical 
field  changes  as  a  function  of  the  dielectric  properties  of 


material,  which  in  turn  changes  the  resonant  frequency. 
Figure  1  depicts  the  overall  configuration  of  the  current 
nondestructive  testing  setup.  The  nondestructive  test 
is  performed  by  evaluating  a  bunch  width  of  reflected 
microwave.  This  is  because  the  increasing  the  whole  of 
cavity  in  cable  insulation  can  be  detected  as  variation 
of  the  bunch  width.  Although  the  method  has  been 
a  conventional  method  for  characterizing  aging  proper¬ 
ties,  there  are  several  practical  disadvantages.  One  is 
that  the  detection  and  characterization  of  variation  of 
bunch  width  are  not  feasible  for  on-site  inspection.  Also 
continuous  surveillance  of  the  variation  is  very  crucial. 
Taking  this  into  account,  our  concern  in  this  paper  is  to 
develop  a  computational  method  for  identifying  dielec¬ 
tric  properties  of  electrical  cables  based  on  a  microwave 
measurement  system. 

The  measurement  of  material  aging  parameters  con¬ 
tinues  to  be  a  very  challenging  problem.  In  our  ap¬ 
proach,  we  formulate  a  forward  problem  arising  in  mi¬ 
cro  wave  NDE.  More  specifically,  it  represents  a  math¬ 
ematical  model  of  a  specific  measurement  using  the 
parameter-to-output  mapping  with  the  appropriate  ad¬ 
missible  class  of  material  parameters,  To  this  end,  the 
mathematical  model  is  described  by  Maxwell’s  equa¬ 
tions  in  2D.  The  numerical  scheme  for  solving  the  for¬ 
ward  problem  is  essentially  the  finite-difference  time- 
domain  method  (FDTD)  in  two  spatial  dimensions. 
Thus,  the  inversion  of  material  parameters  is  imple¬ 
mented  with  the  aid  of  the  forward  problem.  The 
nonlinear  least  square  identification  (OLSI)  is  a  con¬ 
ventional  inverse  methodology  and  there  have  been 
many  efforts  in  a  variety  of  industrial  applications.  For 
OLSI  related  to  electromagnetic  inversion,  we  refer  to 
[3]  [4]  [5]  [6]  [7].  However, it  is  well  known  that  the  prob¬ 
lem  mentioned  above  has  many  solutions  due  to  the  fact 
it  is  ill-posed.  Recently,  interest  has  grown  in  stochastic 
inversion  using  Markov  Chain  Monte  Carlo  (MCMC) 
methods  [8]  [9].  Previous  efforts  based  on  a  similar  ap¬ 
proach  has  been  proposed  [10]  [11].  The  method  has 
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great  advantages  for  the  more  practical  aspects  of  in¬ 
versions  such  as  setting  initial  guesses  and  overcoming 
local  minimums  of  the  inverse  solution.  In  this  paper,  a 
stochastic  inversion  technique  based  on  the  Metropolis- 
Hasting  algorithm  is  applied  to  the  problem  presented 
here. 


Gunn  Oscillator 
(X-band  ca.  9  GHz) 


Attenuator 
/directional  coupler 

Circulator 


Oscilloscope 


FET  Amplifier 


Fig.  1:  Nondestructive  test  using  TRMDA 


2  Forward  Analysis  for  NDE 

2.1  Problem  Domain 

For  the  simplicity  of  discussions,  we  consider  a  two 
dimensional  bounded  domain  x  =  (aq,^)  £  fi  C  R2 
as  an  inspection  area.  Let  fic  be  a  cable  domain  de¬ 
fined  on  fi.  The  cable  domain  consists  of  a  material 
Qd  and  a  metal  core  fi0.  The  material  corresponds  to 
a  polymer  of  which  the  dielectric  constant  er  and  con¬ 
ductivity  <Jd  often  vary  at  different  frequencies.  Fig¬ 
ure  2  depicts  the  two-dimensional  problems  space  con¬ 
sidered  here.  Our  electromagnetic  interrogation  is  then 


Problem  Domain  Q 


Fig.  2:  Schematic  diagram  of  geometry 


represented  by  electric  and  magnetic  fields  in  fi.  Let 


Ds(t,x),  H\(t,x),  H2(t,x)  and  Es(t,x)  be  the  compo¬ 
nent  of  electric  flux  density,  magnetic  and  electric  field 
at  time  t  G  [0,  T]  and  at  location  x  G  fi.  In  order  to  use 
Gaussian  units,  we  also  use  the  normalizing  terms, 


E3  ~  \[^e3 

V  ^0 


where  eo  and  fi o  are  the  permittivity  and  magnetic  per¬ 
meability  of  air,  respectively.  From  Maxwell’s  equa¬ 
tions,  the  electromagnetic  propagation  in  the  transverse 
magnetic  (TM)  mode  is  governed  by 


dD3(t,x) 

1  fdH2(t,x)  dHi(t,x)' 

\ 

dt 

y/eoid0  V  dxx  dx2  , 

) 

+Jg(t,x)  in  [0 ,T]  x  fi 

(1) 

D3(iv,x) 

=  e*(u>)E3(u>,  x) 

(2) 

dHi(t,x) 

 1  dE3(t,  x) 

dt 

,/eo/io  dx2 

in  [0,  T\  x  fi 

(3) 

dH2{t,x) 

 1  dE3(t,  x) 

dt 

VeoMo  dxi 

in  [0,  T]  x  fi 

(4) 

with  zero  initial  states  and  with  the  absorption  bound¬ 
ary  conditions.  The  complex  relative  dielectric  constant 
e*  of  the  polymer  fi^  can  be  represented  by  a  lossy  di¬ 
electric  medium  of  the  form; 

e*  =  er  +  .  d  •  (5) 

jue  o 


On  the  other  hand,  The  value  of  this  parameter  in  fia  = 
fi  —  fic  becomes  e*  =  1.  The  value  of  that  in  fi0  is  also 
given  by 

e*  =  1  +  — 

3ve0 

where  cr0  denotes  the  conductivity  of  the  metal  core. 
The  source  current  Js  is  the  test  signal  in  our  inspection 
process  and  is  given  by 


Js(t,x)  =  5(Qs)gs(t)I(0,t.)(t)  (7) 

where  5  denote  the  array  of  actuator  in  fi  such  as 
fis  p|  fic  =  0.  In  Eq.  (7),  the  test  signal  is  truncated  at 
a  finite  time  ts  by  the  indicator  function  I.  Detections 
can  be  made  through  the  allocation  of  sensor  arrays  fiG 
in  fi  by 

Y(t,x-,q)  =  E3(t,x-,q)  in  0,o  (8) 

where  the  parameter  vector  q  to  be  identified  is  denoted 
by 

q={er  ,cr4.  (9) 

Thus  the  forward  problem  for  our  NDE  is  formulated 
as 

q  ==^  Y (t,  x;  q).  (10) 


2 

Approved  for  public  release;  distribution  unlimited. 


Scattered  Field 


2.2  Numerical  Scheme 

The  finite-difference  time-domain  method  (FDTD)  in 
two  spatial  dimensions  can  be  implemented  to  solve  the 
forward  model  [12].  Putting  Eqs.(l),  (3),  and  (4)  into 
the  finite  difference  scheme  results  in  the  following  dif¬ 
ference  equations: 

-°3+1/2m  ~-pr1/2M 

At 

1  {  H2[i  +  1/2,  j]-H$[i- 1/2,  j\ 

V^oMo  l  Axi 

H?[i,j  +  1/2] 1/2] 

Ax2 

Hi+1[i,j  +  1/2]  —  Hi\i,j  +  1/2] 

At 

_ 1  fB%+1/2[i,j  +  l]-E%+1/2[i,j]\ 

V/eoMo  \  &x2  j 

g2w+1[z,j  +  l/2]-g2"[i,j  +  l/2] 

At 

_ 1  f^3ra+1/2[i  +  l,j]-^3ra+1/2[i,j]f 

V6oMo  [  Axi  J 

In  order  to  prevent  spurious  reflections  from  the  edge  of 
the  problem  space,  the  perfectly  matched  layer  (PML) 
is  adopted  to  two  dimensional  domain  (  See  [13]  for 
more  details  ).  The  principle  of  micro  wave  testing  can 
be  represented  by  the  simulation  of  plane  waves.  Since 
the  propagating  plane  wave  must  not  interact  with  the 
absorbing  boundary  conditions  like  PML,  the  problem 
space  can  be  divided  up  into  two  regions  such  as  total 
fields  and  scattered  field.  Figure  3  illustrates  how  the 
incident  plane  wave  behaves.  To  simulate  a  plane  wave 
interacting  with  a  cable  object,  the  cable  object  accord¬ 
ing  to  its  electromagnetic  properties  must  be  specified. 
The  simulation  of  a  plane  wave  pulse  hitting  a  dielectric 
cylinder  with  er  =  30  and  cr  =  0.3 [m/s]  is  demonstrated 
in  Fig.  4.  After  the  time  step  T  =  50,  the  pulse  is  in¬ 
teracting  with  the  cylinder.  At  time  step  T  =  75  some 
of  it  passes  through  the  cylinder  while  the  other  part  of 
it  goes  around  that.  It  could  be  also  recognized  that, 
at  T  =  100,  the  main  part  of  the  propagation  is  being 
subtracted  out  the  end  of  the  total  field. 


3  Stochastic  Inverse  Analysis 

Our  problem  is  to  identify  dielectric  parameters  re¬ 
lated  to  cable  insulation 


q  =  W,crd} 

using  the  information  obtained  through  a  scattering  ex¬ 
periment  for  the  forward  problem  stated  in  the  previous 
section.  The  observations  consist  of  sampling  data  as¬ 
sociated  with  the  electrical  field.  We  suppose  that  the 
experimental  observations  consist  of  the  values  of  elec¬ 
trical  field  at  the  measurement  point  xv  —  {x^x^}  at 


Fig.  3:  Diagram  of  a  plane  wave  hitting  a  cable  object 


the  sampling  time  interval  At.  Let  Y ™  =  {Yd[n]}^=1 
denote  the  collection  of  the  data  we  seek  to  reconstruct 
and  YiV(q)  =  {Y (nAt,  xv\  q)}^=1  be  the  electrical  field 
arising  from  the  forward  problem  with  dielectric  param¬ 
eter  vector  q.  Then  the  structural  equation  model  is 
formulated  as 

Yd[n\  =  Y(nAt,xp;q)  +r)n,  rjn  ~  .V(0,t2), 

(n  =  1, 2,  •  •  •  N).  (11) 


where  r\n  denotes  the  mutually  independent  sequence 
of  additive  measurement  noise.  With  the  background 
knowledge  of  Bayes  formula,  the  posteriori  density  func¬ 
tion  with  respect  to  the  set  of  unknown  parameter  vec¬ 
tor  q  can  be  simply  written  by 


p{er\a,YN)  oc  p(er)  — =L- 

n=l  ^2?rT 


exp 


Y(nAt,xp;q)  -  Yd[n] 


2t2 


(12) 


p(ad\er,YN)  oc  p{a)  — =L- 

n=l  ^2?rT 


exp 


Y(nAt,xp,q)  -  Yd[n] 


2  r2 


(13) 


where  p(er)  and  p{(Jd)  denote  a  priori  density  functions 
of  er  and  ad.  An  estimation  algorithm  can  be  performed 
by  sampling  procedures  for  the  posteriori  distribution 
given  by  Eqs.  (12)  and  (13)  from  which  sample  paths 
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can  be  drawn  using  Markov  chains.  For  the  practi¬ 
cal  implementation  of  MCMC,  Metropolis-Hasting  al¬ 
gorithm  is  applied  to  the  problem  considered  here. 


(a)  T  =  25 


(b)  T  =  50 


(d)  T  =  100 


4  Experimental  Results 

To  test  the  feasibility  of  the  estimation  approach, 
we  produce  synthetic  data  for  the  observations  = 
{Y[n\ }n=i  by  adding  random  noise  to  the  results  of 
the  simulation  with  a  known  set  of  parameters  qtrue  = 
{e*rne,  crt([ue}.  The  data  to  be  identified  are  generated 
with  the  forward  problem  in  the  previous  section  using 
the  following  parameter  values: 

4™e  =  30.0 

Glpe  =  0.3  [S/m] 

The  carrier  frequency  of  the  plane  wave  is  90 GHz  and 
the  duration  of  sampling  time  is  taken  as  At  =  [s], 
respectively.  The  observation  data  consist  of  512  mea¬ 
surements  of  the  electric  field  taken  at  xp  =  (0.5,  0.5). 
Figure  5  depicts  the  model  data  with  true  parameters. 
A  priori  densities  of  er  and  a  are  preassigned  as  normal 
distributions; 


P(er) 

P(<7d) 


1 

\/27TTr 

1 

427 TTa 


exp 


exp 


V  27* 


where  4  and  cr°  denote  the  nominal  values  of  dielec¬ 
tric  parameters.  The  artificial  noise  term  was  provided 


Fig.  5:  Synthetic  data  with  true  parameter  values  (noise 
free) 


using  a  standard  Gaussian  random  generator  7V(0,  0.1). 
After  10000  iterations  and  discarding  the  first  1000  it¬ 
erations,  the  posterior  means  of  q  were  evaluated.  Typ¬ 
ical  paths  based  on  9000  draws  are  presented  in  Fig.  6. 
The  resulting  marginal  posteriors  appear  in  Fig.  7.  Es¬ 
timation  results  are  summarized  in  Table  1. 


Fig.  4:  Simulation  of  a  plane  wave  pulse  ^ 
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5  Conclusions 

Microwave  analysis  was  considered  for  the  detection 
and  the  characterization  of  cable  insulation.  The  for¬ 
ward  problem  in  2D  was  considered  by  constructing 
deterministic  parameter-to-output  mappings  based  on 
the  frequency  dependent  media  with  material  quanti¬ 
ties  of  cable  insulation.  A  numerical  scheme  was  de¬ 
veloped  using  the  FDTD  method  in  order  to  solve  the 
forward  problem.  Secondly,  stochastic  inverse  method¬ 
ology  was  considered  with  the  background  knowledge 
of  Bayes  formula.  The  full  probability  model  was  spec¬ 
ified  by  the  deterministic  formula  treated  in  the  for¬ 
ward  problem  and  by  the  stochastic  relationship  for 
their  measurement  strategies.  The  posteriori  density 
function  was  provided  with  likelihood  functional  up¬ 
date  based  on  iterative  procedures  of  measurements  and 
with  a  priori  probability  densities  with  respect  to  the 
unknown  quantities.  The  validity  and  feasibility  of  our 
proposed  method  were  demonstrated  through  compu¬ 
tational  experiments  for  appropriate  specific  examples. 
The  MCMC  samplers  using  Metropolis-Hasting  algo¬ 
rithm  was  tested  for  our  inverse  problem. 
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